Universality of Decay out of Superdeformed Bands in the 190 Mass Region 
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Superdeformed nuclei in the 190 mass region exhibit a striking universality in their decay-out 
profiles. We show that this universality can be explained in the two-level model of superdeformed 
decay as related to the strong separation of energy scales: a higher scale related to the nuclear 
interactions, and a lower scale caused by electromagnetic decay. Decay-out can only occur when 
separate conditions in both energy regimes are satisfied, strongly limiting the collective degrees 
of freedom available to the decaying nucleus. Furthermore, we present the results of the two-level 
model for all decays for which sufficient data are known, including statistical extraction of the matrix 
element for tunneling through the potential barrier. 



I. INTRODUCTION 

It is well known that, for a major-to-minor axis ratio 
of about 2, a new set of shell closures and magic num- 
bers occurs in many nuclei. Such superdeformed (SD) 
states are one of the most striking predictions of the sheU 
model [T]. High electric quadrupole moments and small 
centrifugal stretching mark these states as fundamentally 
different from their normally deformed (ND) isomers 
This contrast has stimulated an abundance of experimen- 
tal and theoretical studies, yet several pressing questions 
persist [H, [2, Q. Of these, perhaps the most interesting 
is the mechanism by which SD bands decay. 

After their formation at high angular momentum, typ- 
ically via heavy ion collisions, these nuclei decay to the 
yrast SD rotational band, and then uniformly down that 
band by E2 transitions. The SD bands are observed to 
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FIG. 1: (a) Decay profiles of several SD bands near A ~ 190. 
Note how suddenly each decay-out occurs, (b) The profiles 
of (a), but shifted in angular momentum so that the leftmost 
points, the last point in which the SD band is experimentally 
observed to retain any strength, are aligned. In addition to 
the abruptness of their decay, the profiles are seen to exhibit 
a universal behavior. Both graphs are from Ref. 
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FIG. 2: (Color online) SD branching ratio Fs, calculated in 
the two-level model and showing the onset of universality. 
Universality of the decay profiles arises from the fact that 
interband decay is nearly forbidden until conditions are fa- 
vorable in both energy regimes, at which point SD strength 
vanishes quickly. The sudden transition from intraband to 
interband decay occurs when both Fjv ^ Fs and V > 14, 
where the critical tunneling matrix element Vc is given by 
Eq. (|24p . When both conditions are satisfied, the curves are 
thus nearly identical, giving rise to the observed universality 
of decay profiles. 

retain their strength through many states, even after they 
are no longer yrast, with negligible losses. Then, quite 
suddenly, the SD band loses almost all of its strength 
over just one or two states [see Fig. (Ha)], although the 
nucleus is still well above the SD bandhead. After a series 
of statistical decays through unrelated states, the nuclei 
continue via i?l-dominated decays in an ND rotational 
band 0,3. 

By far the most SD decays have been observed in the 
"classic" 190 mass region. Recently, Wilson and collab- 
orators Q demonstrated a striking feature of SD decay 
in this region: the decay profiles, when corrected for dif- 
fering angular momenta, are nearly identical [Fig. [Ijb)]. 
This universality of decays is to be contrasted with mere 
abruptness, a feature which has been acknowledged for 
some time. Indeed, it was noted as early as Ref. [5] that 
a purely statistical model would be insufficient to ex- 
plain true universality, i.e. strong consistency between 
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different decay profiles. Likewise, a chaos-assisted phe- 
nomenon cannot, of itself, generate universality. 

The purpose of this article is to demonstrate that 
within the two-level model of SD decay-out, the decay 
profile is universal. In this model, the branching ratio is 
completely determined by four parameters: the detuning 
A = En — £s, a tunneling matrix element V, and elec- 
tromagnetically induced broadenings for the SD and ND 
wells, Ts and Tn, respectively. Because nuclear forces 
are very much stronger than electromagnetic, there is a 
strong separation of scales A Fg, r^y. We show that 
the decay-out only occurs when conditions have become 
favorable in both energy regimes. As a consequence, de- 
cay occurs not only very suddenly, but also in a lim- 
ited region of the model's full parameter space. Within 
this subspace, branching ratios are highly insensitive to 
changes in the parameters, resulting in the universality 
observed in experiment. 

Figure [2] illustrates these findings. It depicts the calcu- 
lated value of the in-band branching ratio Fs for various 
ratios Ts/^n- Vc, a simple function of A and Ts/^n 
in the two-level model, sets the scale V must achieve to 
allow decay. Universality is evident from the figure: if 
either Fs/Fat or Vc/V is too large, no decay can occur. 
On the other hand, as both cross critical values, Fs sud- 
denly vanishes. Furthermore, the curves converge in this, 
the decay-allowed limit, resulting in a universal profile in 
good agreement with the experimental results of Fig. [1] 
However, the connection between this theoretical univer- 
sality and the experimentally observed universality still 
needs to be made, because the abscissas of the two plots 
(Figs. 1 and 2) are not the same. This connection will 
require a theoretical understanding of the angular mo- 
mentum dependence of V, which in turn requires a more 
detailed understanding of the barrier between the two 
wells. 

The outline of the paper is as follows. We first briefly 
review the two-level model in Section II. Section III 
presents the model's results, on which Fig. [2] is based. 
Section III further includes a numerical analysis of all 
SD decays for which sufficient data are known. Section 
IV gives our conclusions. 



II. TWO-LEVEL MODEL OF SD DECAY 

Theoretical efforts to describe the SD decay pro- 
cess have centered on a potential function of both nu- 
clear quadrupole deformation and angular momentum. 
Vigezzi and collaborators noted early on that a double 
well in the deformation more accurately models the ex- 
perimental data than any alternative [6|. In this pic- 
ture, the shape of the tunnel barrier and the two wells 
varies as the nucleus sheds angular momentum, and the 
states of the ND and SD wells are broadened by their re- 
spective couplings to the electromagnetic field. Shortly 
thereafter, Khoo et al. Q conclusively demonstrated that 
these electromagnetic widths must be much less than the 
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FIG. 3: Schematic diagram of the two-level model of SD de- 
cay. In each well of the double-well potential, only one level is 
kept. £jv and es are the unperturbed energies of the isolated 
ND and SD states, which are connected by a tunneling matrix 
element V. The two states have electromagnetic decay rates 
rjv/?i and Ts/h, respectively. 

inter-level spacings in each well. The most appropriate 
picture for SD decay is thus found to be two sets of dis- 
crete, slightly broadened states, connected by matrix el- 
ements to tunnel through the barrier. 

The two-level model for SD decay [1] is given by keep- 
ing only one level in each well, the decaying SD level and 
the ND level with nearest energy (see Fig. . Since the 
role of additional ND levels is principally to steal decay 
strength from the first [§, , it is now well established 
that going beyond this level of approximation is not use- 
ful for most heavy- nuclei SD decays. 

A. Green's Function Description of SD Decay 

The Hamiltonian of the two-level model is a sum of 
three terms: H = Hw + Ht + He,. In the basis of the 
two isolated levels, the first term is diagonal: 

where is the energy of state i. H\y generates time 
evolution within each well: its related retarded Green's 
function is Gw{E) ^ {E - Hw + iO+y^ . 
The second term, 

allows tunneling through the barrier. Here we have cho- 
sen the relative phases of the basis states \S) and |A^) 
such that V is positive, without loss of generality. To- 
gether, Hw and Ht form a simple problem common to 
many introductory quantum mechanics texts. 

The remaining term, Hd = Hem + ^c, gives the elec- 
tromagnetic decay of the two levels. Hem is composed 
of the electromagnetic (harmonic oscillator) modes of the 
environment, and He gives their couplings to the nucleus. 
In this case, it is not necessary to treat the particulars 
of these terms; rather we work at the level of the exper- 
imentally determined decay rates Ts/h and Fjv/?i 
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The self-energy due to Hd is 

y__l(Ts 

Dyson's Equation 

G{E) - ([Gw{E)]-^ -Ht~^ 



(3) 



(4) 



gives the Green's function of the fuU system, summing 
the effects of S and Ht to all orders. Treating the physics 
of the two wells {Gw)^ electromagnetic decay (E), and 
the barrier {Ht) on the same footing in this way is essen- 
tial to a complete description of SD decay-out: all three 
play equally important roles in determining experimental 
observables, such as branching ratios. The full Green's 
function of the two-level model is thus 
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[{E-es + irs/2) {E-eN+ *rjv/2) - V^] ' 



E - En + iTN/'2. 
V 



V 

E-es + iTs/2 



(5) 



At t = 0, the nucleus is localized in the SD well by 
virtue of its previous, measurable E2 decay. For later 
times, then, the probability to find the nucleus in the SD 
or ND well is given by Piit) = |Gjs(OP> where S,N 
respectively. Here 

|^G,5(£;)e-^*/^ (6) 

the Fourier transform of Gij{E), is the retarded propa- 
gator from S to i. 

The resulting probabilities are 
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""/'^ [cosh (uJit) ~ cos (l^rt)] (7) 



and 
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F' — huji 



huj, + F' 



ihuJr 



ifwJr 



(8) 



Here, F ee (FAr-|-Fs)/2 and F' = (FAr-Fs)/2, while the 
real and imaginary parts of the complex Rabi frequency 
w are given by 



respectively, where the sign is used for the real part, 
and the "— " for the imaginary. As Eqs. ([7]) and ^ 
demonstrate, LOi is associated with decoherence due to 



coupling with the electromagnetic field, while uj^ is anal- 
ogous to the real Rabi frequency of a closed two-well 
system. 

The branching ratios Fi — (Ti/fi) dtPi[t) are found 
by time-integrating the probabilities. The results are Q 
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Fs = 



FArFV(FAr+F^) 
Fs+FArFi/(FAr+Fi) 



Ts 

Fs + FjvFi/(Fw + Fi)' 



where 



F 



2rv^ 



(lOa) 
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(11) 



Equations pU]) are the expected results for series decay 
out of a two-level problem. In this light, it is clear that 
F^ is simply the net rate for the nucleus, starting in the 
SD well, to tunnel irreversibly through the barrier. 

These results allow us to extract information about 
the potential barrier from experiment. In particular, the 
values determined by a typical SD decay experiment are 
Fg = 1 — F]^ and F5, while Fjv can be estimated by 
applying the cranking model to a Fermi gas density of 
states jj. From Eq. we find 



F^ = Fc 
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Is 
Fjv 



(12) 



B. Determination of V 

To uniquely extract V itself requires A, which in turn 
implies detailed knowledge of the spectrum in the ND 
well. In the absence of this, we consider the entire sta- 
tistical ensemble of two-level models, each characterized 
by a different value of the unknown variable A. 

To proceed, we construct a probability density func- 
tion V{A), which gives the statistical weight each value 
of A has in the ensemble. The simplest ansatz for the 
distribution of energy levels in the ND well is the "struc- 
tureless" Wigner surmise [l3{ : 
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where s is the level spacing in units of its average value 
D]s[. A is the detuning between the SD state and its 
nearest ND neighbor; thus its magnitude must be less 
than half the spacing sD^ between the ND levels just 
above and just below the SD level. Given this spacing, 
therefore, A is drawn from the rectangular probability 
density 
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where 8 is the Heaviside step function. The total prob- 
abihty theorem yields the desired result [§]: 
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(15) 

where erfc(a;) denotes the complementary error function 
of X. 

A probabilistic statement like Eq. (fT5|l obviates the 
need for exact knowledge of A. To arrive at the proba- 
bility density function V(V), one need only perform an 
elementary change of variables: 



V{V) = 2P(A) 



dA{V) 



dV 



(16) 



where the factor 2 results from our choice of phase for V. 
Here, |A(y)| is a function, not the random variable A; 
it is found from Eq. (fTTj) to be 



(17) 



^r^r is the smallest V consistent with 



where Vmin = 
the two-level model. V{V) is thus seen to be [1, [31 
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ViV) = { Dn ri\A{v)\ 




erfc 



: |A(y)| 



otherwise 



(18) 

A probability distribution such as Eq. represents 
the most one can say about V without microscopic knowl- 
edge of the ND well. The mean of 'P{V) is 
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while the standard deviation is 
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(19) 



(20) 



av/{V) « 84%, indicating that VlV) is well peaked 
about (V), and thus this mean provides a good measure 
of the likely value of V. 

The Wigner surmise (fT3|) . provides a reasonable, neu- 
tral guess at the spacings of states in the ND well. It is 
closely related, and may be considered a good approxi- 
mation, to the level distribution of the Gaussian orthog- 
onal ensemble [l^. Nevertheless, it is straightforward to 
reproduce the preceeding analysis, substituting a level- 
spacing density of choice for Eq. p^ . 



III. RESULTS FOR THE 150 AND 190 MASS 
REGIONS 

TableUgives the values of and (V) for all SD decays 
for which the four parameters, F^; Fg; Fjv; and D^, are 



known. In the table, we have further defined the series 
rate to irreversibly leave the SD band: 



/t r It, 



(21) 



It is Tout, directly extractable from experimental results, 
which competes with F5 to determine whether a nucleus 
will decay out of or remain within the SD band. 

The dynamics of SD decay is a consequence of the 
strong separation of energy scales. 



Vs.Tn <^Dn,{V), 



(22) 



exhibited in TableU In the 190 mass region, particularly, 
we find that parameters relating to the potential double- 
well, such as Dm and V, are lO's to lOOO's of electron- 
Volts, while those relating to electromagnetic decay are 
fractions of meV. This is to be expected, since nuclear 
forces are, of course, many orders of magnitude stronger 
than electromagnetic ones. 

Each of these energy scales has a typical rate associated 
with it: oscillations within the two-level system are char- 
acterized by bJr, whereas T jfi gives the typical rate for 
electromagnetic decay in Eqs. (O-®. Since hujr ^ F, 
it is clear that SD decay is primarily a coherent process: 
a nucleus generally undergoes thousands of virtual Rabi 
oscillations during a single decay event. Only if both V 
and A were of order meV or smaller could the decay be 
incoherent, and such "accidental" near-degeneracies are 
masked by the fact that, since F^ « 0, when V IS very 
small, the nucleus cannot leave the SD band. Moreover, 
the probability for A to be of order F is seen from Eq. ([T| 
to be - T/ Dn- 

Equation (jlObp can be rewritten: 



Fs = l- 



1 



1 + (Vc/vy + Fs/Fjv ■ 



where 



(a^ + f') 



Ts/r 



N 



l + Ts/F 



N 



(23) 



(24) 



The critical tunneling matrix element Vc can be extracted 
from experiment via a statistical approach similar to that 
described for V in Section III Bl except that the result 
does not depend on F-^ or Fs- The resulting average 
values {Vc) are tabulated in TableU) 

According to Eq. (|23p , only two dimensionless param- 
eters, Ts/^N and Vc/V, play a role in determining the 
branching ratios; each corresponds to one of the prob- 
lem's two energy scales. It is clear from Table U that 
the first, Fs/Fjv, decreases dramatically over the course 
of each SD band's decay-out. This can be understood 
physically as a relaxation of the nucleus's centrifugal 
stretching, and consequent reduced coupling to the elec- 
tric quadrupole field, as its spin lowers. 

We extracted (V) using the experimental branching 
ratios; thus it would be circular to make use of those 
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TABLE I: Results of the two-level model, for all SD decays for which sufficient data (branching ratios, Ts, Fjv, and Dn) are 
known. / is the nuclear spin quantum number. Note that (14) and Ts/ {Ts + Fjv) do not depend on Fs- The rightmost column 
gives the sources of the experimental inputs and the estimates of Fjv and Dn- 



nucleus(J) 


Fs Ts Tn Dn T^ r„„t (V) (K) Refs. 
(meV) (meV) (eV) (meV) (meV) (eV) (eV) 


"2Hg(12) 
"2Hg(10) 


0.74 0.128 0.613 135. 
0.08 0.050 0.733 89. 


0.049 0.045 8.7 
2.7 0.58 41. 


14.0 0.173 
5.6 0.064 


[4,i5j 


"2pb(16) 
i92pb(14) 
i92pb(12) 

192pb(10) 

i«2pb(8) 


>0.99 0.487 0.192 1,362. 

0.98 0.266 0.201 1,258. 

0.66 0.132 0.200 1,272. 

0.12 0.048 0.188 1,410. 
<0.25 0.016 0.169 1,681. 


<0.0050 <0.0049 <29. 

0.0056 0.0054 34. 

0.10 0.067 170. 

1.9t 0.35 lOOO.t 
>0.067 >0.048 >250. 


288. 0.717 
237. 0.570 
201. 0.398 
160. 0.203 
120. 0.086 


[16, T7] 
[16, TT] 
[16, FT] 
[16, T7] 
[16, jH 


"'Hg(12) 
i3^Hg(10) 


0.58 0.097 4.8 16.3 
<0.09 0.039 4.1 26.2 


0.071 0.070 0.49 
>0.44 >0.40 >2.1 


0.58 0.020 
0.64 0.0094 


[18,i9, 20, 21J 
[18,i9, 20, 21J 


i«4Hg(12) 
"'Hg(lO) 


0.60 0.108 21. 344. 
0.03 0.046 20. 493. 


0.072 0.072 5.0 
1.6 1.5 35. 


6.1 0.0051 
5.9 0.0023 


]22} 
[22] 


"'Hg(lO) 


0.60 0.086 1.345 19. 
<0.05 0.033 1.487 14. 


0.060 0.057 0.97 
>1.1 >0.63 >3.0 


1.2 0.060 
0.52 0.022 


[4, 20] 
[4, 20] 


"'Hg(15) 
''*Hg(13) 
"^Hg(ll) 


0.90 0.230 4.0 26.5 
0.84 0.110 4.5 19.9 
<0.07 0.048 6.4 7.2 


0.026 0.026 0.52 
0.021 0.021 0.34 
>0.71 >0.64 >0.60 


1.5 0.054 
0.77 0.024 
0.15 0.0074 


[20, 21J 
[20, 21] 
[20, 211 


"^Pb(lO) 
i^-^PbCS) 
i94pb(6) 


0.90 0.045 0.08 21,700. 
0.62 0.014 0.50 2,200. 
<0.09 0.003 0.65 1,400. 


0.0053 0.0050 1100. 
0.0087 0.0086 72. 
>0.032 >0.030 >77. 


3300. 0.36 
90. 0.027 
20. 0.005 


[21 , 23 , 24, 25) 
[21,23, 24, 25] 
[21, 23, 24, 25] 


"'*Pb(12) 
"4pb(10) 
i^'*Pb(8) 
i«''Pb(6) 


>0.99 0.125 0.476 236. 

0.90 0.045 0.470 244. 

0.65 0.014 0.445 273. 
<0.04 0.003 0.405 333. 


<0.0013 <0.0013 <2.7 
0.0051 0.0050 6.1 
0.0077 0.0076 8.8 

>0.088 >0.072 >39. 


26.9 0.208 
18. 0.087 
12. 0.031 
7. 0.007 


[17, 2^ 
[17, 2^ 
[17, 2^ 
[17, 21] 


i52Dy(28) 
i^2Dy(26) 


0.60 10.0 17. 220. 
0.19 7.0 17. 194. 


11. 6.7 35. 
140.t 30. 120.t 


33. 0.37 
26. 0.29 


[22] 
[22] 



fCalculated statistically, as explained in the appendix. 



Thus, we explain the universal nature of the decay- 
out profiles as follows: SD decay-out is only allowed 
when suitably low values of both r^/FAr and Vc/V are 
achieved, both of which decrease quickly with decreasing 
spin. The nucleus, therefore, enters the region of allowed 
decay-out very suddenly, moving down the curves from a 
case of very high in-band intensity to one of almost com- 
plete decay-out; high in-band intensity corresponds to 
the long chains of pre-decay SD states observed in exper- 
iment, while the rapid transition into the decay-allowed 
region forms the abrupt decay profiles. However, once 
the system enters the decay-allowed region, the branch- 
ing ratios saturate, and hence are no longer sensitive to 
further changes in the parameters, as Fig. [5] shows. The 
decay profile is consequently universal. 

If ^^^Dy is a representative example, the variation of 
Fg/F^r in the 150 mass region is somewhat slower than in 
the 190 region, although Vc/V still changes dramatically, 
as seen from Fg's approach to its 1^ ^ oo limit. Since 
rapid entry into the decay-allowed region is a necessary 
precondition of universal behavior, we expect that, as 



values in our discussion of Fg^s universality. Instead we 
note that Eqs. (|10b|) and have the limit 

lim Fs = TT^^, (25) 

V^OD 1 s + 1 AT 

and that, in the two-level model, Fs is a monotonically 
decreasing function of V. Values of this limit are given 
in Table HI As we move down each decay chain, it is 
clear that the experimental branching ratios converge to 
these values, and hence we conclude that Vc/V, too, is 
decreasing quickly during each decay chain. This also is 
to be expected: V depends exponentially on the barrier 
height [2^ , which decreases as the ND well drops further 
below the SD one in energy. 

Figure [2] shows Fs as a function of Vc/V and Fs/Fjy. 
As both Fs/Fjv and Vc/V decrease, the SD branching ra- 
tio decreases to an abrupt plateau. Furthermore, within 
the decay-allowed region 

rN>Ts nv>Vc (26) 

all the curves quickly become nearly identical. 
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more data become available for these nuclei, a somewhat 
lesser degree of universality will be observed. 



IV. CONCLUSIONS 

The two-level model has elsewhere been shown 

to be the simplest description of the SD decay-out process 
which still encapsulates the essential physics. It describes 
a two-step decay: first, the nucleus undergoes mainly co- 
herent Rabi oscillations between the SD and ND wells, 
after which it finally decays into one band or the other. 
By using a statistical approach, the two-level model can 
extract as much information as is possible from decay ex- 
periments, including the Hamiltonian matrix element for 
tunneling through the potential barrier, which is of di- 
rect relevance to nuclear structure. Table |T] demonstrates 
the results of this technique for all decays with sufficient 
data, to date. 

Moreover, the most striking property of SD decay-out 
in the 190 mass region, universality of the decay profiles, 
is seen to correspond to universal behavior of the branch- 
ing ratio in the two-level model. The two-level branching 
ratio is completely determined by two dimensionless pa- 
rameters, rs/Fjv and K/F, each of which corresponds 
to one of the two disparate energy scales of the problem. 
Decay-out only occurs when both of these parameters 
are decreasing rapidly. Thus, in the sector of parameter 
space for which decays are allowed, the branching ratios 
saturate and are insensitive to variations of the param- 
eters. Assuming that the relationship between nuclear 
spin and barrier shape is reasonably consistent from nu- 
cleus to nucleus, the resulting decay profiles are neces- 
sarily similar. 



In only two decays of Table [H I92pb(l0) and i52Dy(26), 
is this condition violated. While it is possible that this 
is due to a breakdown of the two-level approximation in 
these cases, in the absence of a physical argument for 
the near degeneracy of two or more ND levels, it is more 
likely that one or more of the input parameters is poorly 
known. T^, in particular, is difficult to estimate, with 
uncertainty crrjv ~ T^r. 

Thus, we estimate statistically for these two decays, 
assuming the true T n differs from the estimated value 
by a "cut" normal distribution: 



ViTN) 



Fat >r„ 



(28) 



otherwise 
where the constant of renormalization due to the con- 
straint is 



A=2i erfc 



1 



V2 



N 



(29) 



Assuming that the two-level approximation is valid, the 
probability density function of follows: 



P(Fi)-P(Fjv) 



dT 



N 



= ( TT ) ^(r^)' (30) 



where F^ is the function of Fat given by Eq. (fT2|) . For 
i92pb(10) and i52Dy(26), TableUgives the median of this 
distribution as the typical value of F^ , from which {V) is 
found. 
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Equation (jl2p places a limit on the experimentally- 
determined quantities. Positivity of F^ requires that 

Eat > Tout (27) 
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